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Extending the concept of parton densities onto nonforward matrix elements (p' j 0(0, z) \ p) 
of quark and gluon light-cone operators, one can use two types of nonperturbative functions: 
double distributions (DDs) f(x,a;t),F(x,y;t) and skewed (off&nonforward) parton distributions 
^/-^ . (SPDs) H(x, £; t), J-((X, t). We treat DDs as primary objects producing SPDs after integration. 

^ ' We emphasize the role of DDs in understanding interplay between X (x) and £ (£) dependences 

\Q , of SPDs. In particular, the use of DDs is crucial to secure the polynomiality condition: iVth 

' moments of SPDs are iVth degree polynomials in the relevant skewedness parameter £ or £. We 

propose simple ansatze for DDs having correct spectral and symmetry properties and derive model 
expressions for SPDs satisfying all known constraints. Finally, we argue that for small skewedness, 
one can obtain SPDs from the usual parton densities by averaging the latter with an appropriate 
weight over the region [X — £, X] (or [x — £, x + £]). 



PACS number(s): 12.38.Bx, 13.60.Fz, 13.60.Le 
<D , I. INTRODUCTION 
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Nonforward matrix elements (p — r | 0(0, z) \ p) \ z 2 =0 of quark and gluon light-cone operators which appear in 
k>( \ applications of perturbative QCD to deeply virtual Compton scattering (DVCS) and hard exclusive electroproduc- 
tion processes Q-5lcan be parametrized by two basic types of nonperturbative functions. The double distributions 
(DDs) F(x,y;t) [2]3j specify the Sudakov light-cone "plus" fractions xp + and yr + of the initial hadron momentum 
p and the momentum transfer r carried by the initial parton. Treating the proportionality coefficient £ as an in- 
dependent parameter one can introduce an alternative description in terms of the nonforward parton distributions 
(NFPDs) J-((X; t) with X = x + yC, being the total fraction of the initial hadron momentum taken by the initial 
parton. The shape of NFPDs explicitly depends on the parameter £ characterizing the skewedness of the relevant 
nonforward matrix element. This paramctrization of nonforward matrix elements by >Ff(X;t) is similar to that 
proposed by X. Ji ]jj who introduced off-forward parton distributions (OFPDs) H(x, £;t) in which the parton 
momenta and the skewedness parameter £ = r + / P + are measured in units of the average hadron momentum 
P = (p + p')/2. There are one-to-one relations between OFPDs and NFPDs M, so it is convenient to treat them 
as particular forms of skewed parton distributions (SPDs). 

In our approach, DDs are primary objects producing SPDs after an appropriate integration. Our main goal in 
this letter is to show that using the formalism of DDs (in particular, their support and symmetry properties) one 
can easily establish important features of SPDs such as nonanalyticity at border points X — £, and x = ±£, 
polynomiality of their X N and x N moments in skewedness parameters £ and £, etc. We also discuss simple models 
for DDs which result in realistic models for SPDs. 



*Also Laboratory of Theoretical Physics, JINR, Dubna, Russian Federation 



II. DOUBLE DISTRIBUTIONS AND THEIR SYMMETRIES 



In the pQCD factorization treatment of hard electroproduction processes, the nonperturbative information is 
accumulated in the nonforward matrix elements (p — r | O(0, z) \ p) of light cone operators 0(0, z). For z 1 = the 
matrix elements depend on the relative coordinate z through two Lorentz invariant variables (pz) and (rz). In 
the forward case, when r = 0, one obtains the usual quark helicity-averaged densities by Fourier transforming the 
relevant matrix element with respect to (pz) 

(p, s' | M0)zm *\ A)1> a {z) \p,s)\ Z 2 =0 = u(p, s')zu(p, s) J (e-^f a {x) - / a (x)) dx , (1) 

where E(0, z; A) is the gauge link, u(p', s'), u(jp, s) are the Dirac spinors and we use the notation 7 Q z Q = z. In the 
nonforward case, we can use the double Fourier representation with respect to both (pz) and (rz): 

(p', s' | # a (Q)zE{0, z- A)M*) I P, 5} I *==o (2) 
= u(p', s')zu(p, s) [ dy [ e -i*(pz)-iv(rz) y . ^ 0(o< x + y < 1) dx + il K" - term , 



where the "if" -term stands for the hadron helicity-flip part yJJ^]. For any Feynman diagram, the spectral constraints 
— 1<x<1,0<j/<1,0<x + 3/<1 were proved in the a-representation || using the approach of Ref. ||. The 
support area for the double distribution F a (x,y;t) is shown on Fig. [Ha. 
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FIG. 1. a) Support region and symmetry line y = x/2 for j/-DDs F(x,y\t); b) support region for a-DDs f(x,a). 

Comparing Eq. (|lj) with the r = limit of the DD definition @ gives the "reduction formulas" relating the 
double distribution F a (x,y;t = 0) to the quark and antiquark parton densities 



F a (x, y; t = 0)| x>0 dy = f a (x) 



F a (x,y;t = 0)\ x<0 dy = -fa{-x) ■ 



(3) 



Hence, the positive-x and negative-x components of the double distribution F a (x, y; t) can be treated as nonforward 
generalizations of quark and antiquark densities, respectively. If we define the "untilded" DDs by 



F a (x,y;t) = F a (x,y;t)\ x>0 ; Fa(x,y;t) = -F a (-x, 1 - y,t)\ x<0 , 

then x is always positive and the reduction formulas have the same form 

-X 

F at a(x,y;t = 0)1^0 <*/ = fa,a( x ) 



(4) 



(5) 

in both cases. The new antiquark distributions also "live" on the triangle < x, y < 1, < x + y < 1. Taking z 
in the lightcone "minus" direction, we arrive at the parton interpretation of functions F aj a(x,y;t) as probability 
amplitudes for an outgoing parton to carry the fractions xp + and yr + of the external momenta r and p. The 
double distributions F(x 1 y;t) are universal functions describing the flux of p + and r + independently of the ratio 
r + /p + . Note, that extraction of two separate components F a (x, y; t) and F a (x, y; t) from the quark DD F a (x, y; t) 
as its positive-x and negative-x parts is unambiguous. 



2 



xp+yr 



xp-(l-y) r 




p-r P 





p-r p 



yri Ul-y)r 



I 



xP+(l+a) r/2 



X 



p-r P+r/2 




xP-(l-a)r/2 



P-r/2 



") b) c) d) 

FIG. 2. a) Parton picture in terms of y-DDs; b, c) -FW-type contributions; d) parton picture in terms of a-DDs. 

In principle, we cannot exclude the third possibility that the functions F(x,y;t) have singular terms at x = 
proportional to 5(x) or its derivative(s). Such terms have no projection onto the usual parton densities. We 
will denote them by FM(x,y;t) — they may be interpreted as coming from the i-channel meson-exchange type 
contributions (see Fig.^jb). In this case, the partons just share the plus component of the momentum transfer r: 
information about the magnitude of the initial hadron momentum is lost if the exchanged particle can be described 
by a pole propagator ~ l/(t — m 2 M ). Hence, the meson-exchange contributions to a double distribution may look 
like 



F+(x,y-,t)~6(x)^&- or F~ (x, y; t) ~ 6'(x) ^ {y) 



M 



t 



M 



f 



etc. 



(6) 



where f^iv) are ^ ne functions related to the distribution amplitudes of the relevant mesons M ± . The two examples 
above correspond to x-even and x-odd parts of the double distribution F(x,y;t). The singular terms can also be 
produced by diagrams containing a quartic pion vertex (Fig.^c) Q] . 

To make the description more symmetric with respect to the initial and final hadron momenta, we can treat 
nonforward matrix elements as functions of (Pz) and (rz), where P = (p + p')/2 is the average hadron momentum. 
The relevant double distributions f a {x, a ; t) [which we will call a-DDs to distinguish them from y-DDs F(x, y; t)] 
are defined by 



(p'\$ a (-z/2)zMz/2)\p) = u(p')zu(p) 



1 rl-\x\ 

dx / e" 

-1 J-l+\x\ 



■«(^*)-<o(r*)/a f a (x,a\t)da + "fc" -terms. 



(7) 



The support area for f a (x,a;t) is shown in Fig.|l|b. Again, the usual forward densities f a (x) and fg.(x) are given 
by integrating f a (x, a; t = 0) over vertical lines x = const for x > and x < 0, respectively. Hence, we can split 
f a (x,a; t) into three components 



f a (x, a;t) = f a (x, a ; t) 9{x > 0) - / a (-x, -a ; t) 9{x < 0) + f M {x, a;t), 



(8) 



where Jm{x 7 a ; t) is a singular term with support at x = only^. Due to hermiticity and time-reversal invariance 
properties of nonforward matrix elements, the a-DDs are even functions of a: 

f a (x,a;t) = f a (x,-a;t) . 

For our original y-DDs F a a (x,y;t), this corresponds to symmetry with respect to the interchange y «-> 1 — x — y 
established in Ref. In particular, the functions <p^{y) for singular contributions Fj^(x, y;t) are symmetric 
^iiv) = ^Af (1 ~ V) both for x-even and x-odd parts. The a-quark contribution 

O s a (-z/2,z/2) = l -[i> a {-z/2)zE{-z/2,z/2; A)^ a {z/2) - {z -» -z}] 

into the flavor-singlet operator can be parametrized either by j/-DDs F^(x, y; t) or by cv-DDs f^(x, a ; t) 



^As argued by M. Polyakov and C. Weiss in the case of pion distributions it makes sense to write the (Pz)-independent 
terms as a separate integral over a single variable y rather than to include them into a singular part of DDs. 
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(p',s'\ O s a (-z/2,z/2) \p,s)\ z 2 =0 

= u(p',s')zu(p,s) [ dx [ \ C e -te(f«)-i(w-i/2)(r») _ jx(pz)+i{ y -i/2){rz)\ p s fa y .Q d y + "jqf " -term 
Jo Jq 2 V / 

= u( P ',s')zu(p,s) [ dx f ' e-" (Pz) - M(rz)/2 /,f(a;,a; t) da + "fcf'-term. (9) 
J-i ./-i+|x| 

In the second line here we have used the fact that positive- x and negative- a; parts in this case are described by the 
same untilded function 

F a{ x -,Vit)\x^0 = F a (x,y;t) + F a (x,y\t). 
The a-DDs /,f (x, a ; t) are even functions of a and, according to Eq. (H), odd functions of x: 

f^(x,a;t) = {f a (\x\,\a\;t) + f- a (\x\,\a\;t)}sign(x) + fi(x,a;t) . (10) 
Finally, the valence quark functions f& (x, a > t) related to the operators 

OY(-z/2, z/2) = ^$ a (-z/2)zE(-z/2, z/2; A)^ a (z/2) + {z - -z}] 

are even functions of both a and x: 

/^(x, a; *) = / Q (|x|, \a\-t) - / s (|x|, |a|; t) + f^(x, a; t) . (11) 

III. PARTON INTERPRETATION AND MODELS FOR DOUBLE DISTRIBUTIONS 

The structure of the integral (^|) relating double distributions with the usual ones has a simple graphic illustration 
(see Fig.||): integrating DDs over a line orthogonal to the x axis, we get f(X). 

The reduction formulas and interpretation of the x-variable as the fraction of the p (or P) momentum suggest 
that the profile of F(x,y) (or f(x,a)) in x-direction is basically determined by the shape of /(x). On the other 
hand, the profile in y (or a) direction characterizes the spread of the parton momentum induced by the momentum 
transfer r. In particular, since the a-DDs /(x, a) are even functions of a, it make sense to write 

f(x,a) = h(x,a)f(x), (12) 

where /i(x, a) is an even function of a normalized by 

l-x 

h(x, a) da = 1. (13) 

-l+a; 

We may expect that the a-profile of h(x, a) is similar to that of a symmetric distribution amplitude (DA) (p(a). 
Since |a| < x, to get a more complete analogy with DAs, it makes sense to rescale a as a = xf3 introducing 
the variable (3 with x- independent limits: — 1 < j3 < 1. The simplest model is to assume that the profile in 
the /3-dircction is a universal function g{(3) for all x. Possible simple choices for g(f3) may be 8(0) (no spread 
in /3-direction) , |(1 — (3 2 ) (characteristic shape for asymptotic limit of nonsinglet quark distribution amplitudes), 
y|(1 — /3 2 ) 2 (asymptotic shape of gluon distribution amplitudes), etc. In the variables x, a, this gives 

hM M = 5(a) , /»(*>(*,«) = , ^H X ,a) = ^ { -^-^ . (14) 

It is straightforward to generalize these models onto the "tilded" DDs f(x,a) with x ranging between —1 and 1: 
/(x, a) should be even in x for the gluon and nonsinglet quark distributions and odd in x for the singlet quark 
case. Furthermore, one can construct ansatze for functions f(x,a;t) involving nonzero t values, e.g., the model 

fi(x,a;t) = fo(x,«)/j(x)exp j ^^_" 2 ^ j (15) 

with h(x, a) — 5(a) and experimental valence densities /^(x) was used in ref. Q to describe the Fi(t) form factor 
and wide-angle Compton scattering. 
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FIG. 3. Integration lines for integrals relating SPDs and DDs. 



IV. RELATIONS BETWEEN DOUBLE AND SKEWED DISTRIBUTIONS 

An important parameter for nonforward matrix elements is the coefficient of proportionality £ = r + /p + (or 
£ = r + 1 P + ) between the plus components of the momentum transfer and initial (or average) hadron momentum. 
It specifies the skewedness of the matrix elements. The characteristic feature implied by representations for double 
distributions [see, e.g., Eq.(|2|)] is the absence of the £-dependence hi the DDs F(x, y) and ^-dependence in f(x, a). 
An alternative way to parametrize nonforward matrix elements of light-cone operators is to use £ (or £) and the 
total momentum fractions X = x + yC (or x = x + £a) as independent variables. If we require that the light- 
cone plus components of both the momentum transfer r and the final hadron momentum p — r are positive, then 
< £ < 1 and < £ < 1. Using the spectral property 0<x + y<loi double distributions, we obtain that the 
NFPD variable X satisfies the "parton" constraint < X < 1. Integrating each particular double distribution over 
y gives the nonforward parton distributions || 

/>1 nX—X 

F ( (X)= dx 6(x + ty-X)F i (x,y)dy 

Jo Jo 

/■X/C fX/Q 

= 9(X>C) F i (X-yC,y)dy + e(X<0 F i (X-yC,y)dy, (16) 

Jo Jo 

where £ = 1 — C- The two components of NFPDs correspond to positive (X > £) and negative (X < £) values of the 
fraction X' = X — ( associated with the "returning" parton. As explained in refs. the second component can 
be interpreted as the probability amplitude for the initial hadron with momentum p to split into the final hadron 
with momentum (1 — Q)p and a two-parton state with total momentum r — (p shared by the partons in fractions 
Yr and (1 — Y)r, where Y — X/C,. 

The relation between NFPDs and DDs can be illustrated on the "DD-life triangle" defined by < x, y, x + y < 1 
(see Fig.||a). Specifically, to get Tq{X), one should integrate F(x,y) over y along a straight line x — X — Cv- 
Fixing some value of £, one deals with a set of parallel lines intersecting the cc-axis at x — X. The upper limit of 
the y-integration is determined by intersection of this line either with the line x + y = 1 (this happens if X > C) or 
with the y-axis (if X < £). The line corresponding to X = £ separates the triangle into two parts generating the 
two components of the nonforward parton distribution. 

In a similar way, we can write the relation between OFPDs H (x, £; t) Q and the a-DDs f(x, a; t) 

H(x,^;t) — / dx 6(x + £a — x) f(x,a;t) da . (17) 

J-l J—l+\x\ 

We use here the tilded notation H(x,£,;t) to emphasize that OFPDs as defined by X. Ji (!]] correspond to 
parametrization of the nonforward matrix element by a Fourier integral with a single common exponential. Note 
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that Eq. ( 117] ) allows to construct H(x,£;t) both for positive and negative values of £. Since DDs f(x,a;t) are 
even functions of a, the OFPDs H(x,^;t) are even functions of £: 

H(x,t;t)=H(x,-{;t). 

This result was originally obtained by X. Ji [Tcfl who directly used hermiticity and time reversal invariance properties 
in his definition of OFPDs. 

The delta-function in Eq.(p7|) specifies the line of integration in the {x,a} plane (see Fig.||). For definiteness, 
we will assume below that £ is positive. The integration line x = x — £a consists of two parts corresponding to 
positive and negative values of x. In the case of quarks with flavor a, substituting f a (x, a) by f a (x, a) or fa(x, a), 
respectively (see Eq.(^)), we get OFPD H a (x,^t) as the sum of three components 

H a (x, £ t) = H a (x, £ t) <x<l)- H s (-x, £; t) 6{-l <£<£) + H M (x, £; t) 0(S <x<0, (18) 
where Hm(x,£', t) comes from integration of the singular term Jm(x — £a, a) over x/£ — e < a < x/£ + e and 



H a , s {x, C; t) = 6>(e < x < 1) / i _ f aA (x - £a, a) da + 0(-f f a ,a{x - fa, a) da . (19) 

The OFPD if a (i,^; t) is in a one-to-one correspondence with the "tilded" NFPD T"(X) introduced in our paper 
H . It parametrizes the nonforward matrix element of the quark operator ?/> a . . . tp a through a Fourier integral with 
a single common exponential. The support of Fg(X) is— 1 + C < X < 1 and by 

£%(x) = T£(x) 9(o < x < i) - t^(c - x) 0(-i + c < x < c) + ^ C M (X) e»(o < x < c) (20) 



it is related to the untilded components given by Eq. (jlq) . In the middle region < X < £, the components 
J^' a (X) appear only through the difference fg(X) — J 7 "^— X). In a recent paper JllJ , Golec-Biernat and Martin 

argued that the decomposition of Pf(X) in the middle region into ^Fq(X) and — X) parts made in Ref. j^] 

amounts to "doubling the quark degrees of freedom" ^ Compared to Ref. ||, we have an extra function Tq 1 (X) 

in Eq.(|20|), so one may question now whether it make sense to represent THX) as a sum of three functions in 
the < X < £ region. Of course, if there were only one value of £ in the nature, one would never get an idea 
about how much of J^(X) should be attributed to J^(X),J^(X) or Fg*(X). The crucial missing element is the 

interplay between ( and X dependences. We stress that our decomposition of !F^(X) is based on the splitting of 
the underlying y-DDs F a (x,y) into positive-x, negative- x and zero- a; parts. The DDs contain information about 
NFPDs for all possible £'s and X's, and that is why the DDs produce an unambiguous decomposition: DDs 
"know" not only what is the shape of T^(X) for a particular £, but also how this shape would change if one 
would take another £. The simplest illustration of interplay between X and £ dependences is provided by NFPDs 
F¥ (X) = 9(0 < X/( < l)ip(X/Q)/\C,\ corresponding to singular parts of DDs. Clearly, knowing (X) at some 
C = Co, w e can obtain its shape for any other £ by rescaling. To demonstrate that NFPDs with such a behavior 
can be obtained only from singular DDs, we write a formal inversion of the basic relation ( |l6|) 

/oo poo 
dX A{X-x-(y}F c (X)dC , (21) 
-oo J — oc 

where the (mathematical) distribution A(z) is defined by 

AM = ^/>|e^. (22) 
Taking J 7 ^ 1 (X) = 9(0 < X/( < l)ip(X/()/\(\ and using the following property of the A-function 



'"They also proposed to split our function T^(X) into overlapping < X < 1 and — 1 + C < X < (, parts to introduce 
"off-diagonal" "quark" f"(X) = 6(0 < X < l)f"(X) and "antiquark" f%(£-X) = -0(0 <X< l)T"(C~X) distributions 
both of which include the same middle part of ^(X). 
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/oo 
A[a-(b]d( = 8(a)8(b) ) 
-OO 



(23) 



we obtain from Eq. ( |2l| ) that F M (x,y) — 5(x)ip(y). 

Thus, information contained in SPDs originates from two physically different sources: meson-exchange type 
contributions J-^ 1 (X) coming from the singular x = parts of DDs and the functions ^(X), T^(X) obtained 
by scanning the x ^ parts of DDs F a (x,y), F a (x,y). The support of exchange contributions is restricted to 
< X < (. Up to rescaling, the function J 7 ^ 1 (X) has the same shape for all (. For any nonvanishing X, 
these exchange terms become invisible in the forward limit £ — > 0. On the other hand, interplay between X and 
£ dependences of the functions !F^(X), T"(X) is quite nontrivial and their support in general covers the whole 
< X < 1 region for all £ including the forward limit £ = in which they convert into the usual (forward) densities 
f a (x), f a (x). The latter are rather well known from inclusive measurements. Hence, information contained in 
f a (x), f a (x) can be used to restrict the models for J^(X), JF?(A). Note that the functions F a (x,y) and F a (x,y) 
are independent as are their £-sensitive scans FHX) and Tq(X). Instead of F a (x,y) and F a (x,y), one can use 
as independent functions their sum F a (x,y) + F a {x,y) which contributes to the quark singlet functions and the 
difference F a (x, y)—F a (x, y) which appears in the valence functions. Extending the DDs onto the whole — 1 < x < 1 
segment does not require extra dynamical information: one should only take into account that the singlet term 
f^(x,a) must be odd in x (see Eq.(|To|)) while the valence term fY(%, a) must be even in x (see Eq.([n])). As a 
result, the singlet contribution T^' a (X) is an odd function of X — £/2 while the valence one J-^' a {X) is an even 
function of X — £/2. 

In our approach, DDs are the starting point while SPDs are derived from them by integration. However, even if 
one starts directly with SPDs, the latter possess a property which forces the use of double distributions. According 
to Eq.(|l6|), the X N moment of T^X) must be a polynomial in £ of a degree not larger than N . A similar statement 
holds for off-forward distributions H(x,£,;t): their x N moments are Ath order polynomials of £. As explained by 
X. Ji jK^, this restriction on the interplay between x and £ dependences of H(x,£;t) follows from a simple fact 
that the Lorentz indices /ii . . . ixn of the nonforward matrix elements of a local operator O^i - mn can j-, e carr i ec l 
either by P^* or by . As a result, 

N /N\ N /N\ 

(P r/2\m(d + ) N m\P + r/2) =£ J (P+) N - k (r+) k A m = (P+) N £ j£, k A Nk , (24) 

k=0 ^ ' k=0 ^ ' 

where (^f) = Nl/(N — k)\k\ is the combinatorial coefficient. Our derivation ( |l7| ) of OFPDs from a-DDs automati- 
cally satisfies the polynomiality condition (|24l), since 



f(x 1 a)x 1 ^ k a k da . (25) 

-l+|x| 



Hence, the coefficients A^k m ( p4[ ) are given by double moments of a-DDs. This means that modeling SPDs 
one cannot choose the coefficients Amu arbitrarily: symmetry and support properties of DDs dictate a nontrivial 
interplay between N and k dependences of Ajvfc's. After this observation, the use of DDs is an unavoidable step in 
building consistent parametrizations of SPDs. 

The formalism of DDs also allows one to easily establish some important properies of skewed distributions. Notice 
that due to the cusp at the upper corner of the DD-life triangle, the length of the integration line nonanalytically 
depends on X for X = £. Hence, unless a double distribution identically vanishes in a finite region around the upper 
corner of the DD support triangle, the A-dependence of the relevant nonforward distribution must be nonanalytic 
at the border point X = (. Furthermore, the length of the integration line vanishes when X — > 0. As a result, the 
components !F^' a (X) vanish at X — if the relevant double distribution F a,a (x 1 y) is not too singular for small 

x. The combined contribution of JF?(A) and T"(C, — X) into the total function !F®(X) in this case is continuous 
at the nonanalyticity points X — and X = C- As emphasized in Ref. 0, because of the 1/A and l/(X — £) 
factors (l/(i ± £) factors if OFPD formalism is used) contained in hard amplitudes, this property is crucial for 
pQCD factorization in DVCS and other hard electroproduction processes. Note, that there is also the exchange 
contribution J-^{X). If it comes from a S(x)ip(y) type term and <p(y) vanishes at the end-points y = 0, 1 the 
Fq 1 (X) part of NFPD vanishes at X — and X = £. The total function !F^(X) is then continuous at these 

nonanalyticity points (OFPDs H(i, £; t) in this case are continuous at x = ±£). In the quark singlet case, the DDs 
should be odd in x, hence the singular term involves 6'(x)(p(y) (or even higher odd derivatives of S(x)). One can 
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get a continuous SPD in this case only if ip'(y) vanishes at the end points. Such a restriction might be too strong 
to be satisfied in all cases. In particular, an essentially discontinuous behavior of singlet quark OFPDs for x = ± £ 
was obtained in a nonperturbative (chiral soliton) model fl2| . 



V. MODELS FOR SKEWED DISTRIBUTIONS 



The properties discussed above can be illustrated by SPDs constructed using simple models of DDs specified in 
Section III. In particular, for the model F^°'(x,y) = S(y — x/2)/(x) (equivalent to f(°'(x,a) = S(a)f(x) ), we get 

i.e., NFPDs for non-zero £ are obtained from the forward distribution f(X) = J-^ = o(X) by shift and rescaling. 
This is an example of a peculiar case of a DD with an empty upper corner: it gives NFPDs with no explicit 
nonanalyticity at X = £. "As a compensation", !F^°\X) vanishes not only for X = 0, but on the finite segment 
< X < C/2. Using the relations J| 

H(x,z-,t)\ i> t = (i-<;/2)r c (x-,t)\ x>c ; x = X x Zal > ^ = ^c (27) 

between our nonforward distributions in the X > £ region and Ji's off-forward parton distributions H(x,£;t) jlj 
in the x > £ region, one can see that the narrow F^(x,y) ansatz gives the simplest model H^(x,£;t = 0) = 
f(x) in which OFPDs at t = have no ^-dependence. This result can be obtained directly by using the model 
f(°)(x,a) — <5(a) f{x) for the a-DDs. Another example is the model fll|[l4| in which NFPDs do not depend on £, 
i.e., Tq{X) = f(X). Using the inversion formula ( |2l| ) and Eq. (p3|), we obtain F(x,y) — 5(y)f(x), i.e., the support 
of this DD is on the y-axis only, which violates the mandatory y <-> 1 — x — y symmetry. Unlike the ^-independent 
ansatz for OFPDs, the ^-independent ansatz for NFPDs is forbidden. 

In case of two other models, simple analytic results can be obtained only for some explicit forms of /(x). For 
the "valence quark" -oriented ansatz /W(x, a), the following choice 

f w (x)=Ax- a (l-x) 3 (28) 

is both close to phenomenological quark distributions and produces a simple expression for the double distribution 
since the denominator (1 — a;) 3 factor in Eq. ( |l4| ) is canceled. As a result, the integral in Eq. (|l7]) is easily performed 
and we get 

H^(x, 0\\s\>i = A { [(2 - a)£(l - x)(x\- a + x\- a ) + (£ 2 - i)(x\- a - *r a )] 0(x) + (x -> -x)} (29) 
for \x\ >£, wherei = 3yir(l-a)/2r(4-a), and 

H( lv Hx,0\m<z = A {xl- a [(2 - a)£(l - x) + (£ 2 - x)] + (x -> -$)} (30) 

in the middle — £ < x < £ region. We use here the notation x\ = (x + £)/(l + £) and X2 — (x — £)/(l — £) [^o| . 
As expected, these expressions are explicitly non-analytic for x — ±£. It is interesting to note that in a particular 
case a = 0, the x > £ part of OFPD has the same x-dependence as its forward limit, differing from it by an overall 
^-dependent factor only: 

H^)(x, 0\a=o = A 6(\x\ >0+A 0(\x\ < £) . (31) 

To extend this expression onto negative values of £, one should substitute £ by |£|. One can check, however, that 
no odd powers of |£| would appear in the x N moments of H^ iv \x,^). 

For the singlet quark distribution, the a-DDs f s (x,a) should be odd functions of x. Still, we can use the 
model like (§8]) for the x > part, but take f^ is \x,a)\ x ^ = f^(\x\, a) sign(x). Note, that the integral (|l7|) 
producing H s (x, £) in the |x| < £ region would diverge for a —* x/£ if a > 1, which is the usual case with standard 
paramctrizations of singlet quark distributions at sufficiently large Q 2 . However, due to the antisymmetry of 
f s (x,a) wrt x — > —x and its symmetry wrt a — > —a, the singularity at a = i/£ can be integrated using the 
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principal value prescription which in this case produces the x — > — x antisymmetric version of Eqs.(29) and (|3£ 
As far as a < 2, the resulting functions are finite for all x and continuous at x — ±£. For a — 0, the middle part 
reduces to 

H (1S) (iU)l| £ |<^o = A * %^+\q/ ■ ( 32 ) 

Evidently, the use of the principal value prescription is equivalent to imposing a subtraction procedure for the 
divergent second integral in Eq. ( |19|) defining the untilded functions H a ' a (x,t;). 

In general case, to study the deviation of skewed distributions from their forward counterparts for small £ (or 
£), let us consider the integral producing the x > £ part of H(x, £) [see Eq.(|l9|)1 and expand it in powers of £: 



! /■(!-*) d *f(x,a) 2 a f df(x.a) df(x,a) 
' a da + (1 — x) [ 2- 



2 ./_(i_x) dx 2 \ da dx 



(33) 



where f(x) is the forward distribution. For small £, the corrections are formally 0(£ 2 ), i.e., they look very small. 
However, if f(x, a) has a singular behavior like x~ a , then 

8x 2 x 2 W'^' 

and the relative suppression of the first correction is 0(^ 2 /x 2 ). The corrections are tiny for all x except for the 
region x ~ ^ where the correction has no parametric smallness. Nevertheless, even in this region it is suppressed 
numerically, because the a 2 moment is rather small for a distribution concentrated in the small-a region. It is easy 
to write expicitly all the terms which are not suppressed in the x ~ £ — > limit 

where the ellipses denote the terms vanishing in this limit. Due to strong numerical suppression of higher terms, 
the series converges rather fast. For small x, we can neglect the x-dependence of the profile function h(x, a) in Eq. 
( |l2| ) and take the model f(x,a) — f(x)p(a) with p(a) being a symmetric weight function on —1 < a < 1 whose 
integral over a equals 1. In the region where both x and £ are small, we can approximate Eq. (|l7]) by 

H(x; = J f{x - ta)p(a) da + . . . , (35) 

i.e., the OFPD H(x;£) is obtained in this case by averaging the usual (forward) parton density f(x) (extended 
onto — 1 < x < 1) over the region x — I; < x < x + £ with the weight p(a). In terms of NFPDs, the relation is 

F C (X) = J ^ f(X - C(l + a)/2)p(a) da + . . . , (36) 

i.e., the average is taken over the region X — ( < x < X . 

The imaginary part of hard exclusive meson electroproduction amplitude is determined by the skewed dis- 
tributions at the border point. For this reason, the magnitude of J^(C) [or £)], and its relation to the 
forward densities f(x) has a practical interest. Assuming the infinitely narrow weight p(a) = 5(a), we have 
^{X) = f(X — C/2) + . . . and H(x,£) — f(x). Hence, both ^(C) and H(£,£) are given by f(xBj/2) because 
C = x Bj and £ = xsj/2 + . . .. Since the argument of f(x) is twice smaller than in deep inelastic scattering, this 
results in an enhancement factor. In particular, if f(x) ~ x~ a for small x, the ratio ^c(C)//(C) 1S 2°- The use of a 
wider weight function p(a) produces further enhancement. For example, taking p(a) = f (1 — a 2 ) and f(x) ~ x~ a 
we get 

^ = 1 (37) 

/(C) (l-a/2)(l-a/3) 1 > 

which is larger than 2 a for < a < 2. Due to evolution, the effective parameter a is an increasing function of Q 2 . 
As a result, the above ratio slowly increases with Q 2 . 
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Finally, I want to point out that possible profiles of f(x,a) in the a-direction are restricted by inequalities (see 
|jl4] , fL0 15 lq|) relating skewed and forward distributions. For quark OFPDs, I obtained 15 



v^ 



If one uses the infinitely narrow model f(°'(x, a) = f(x)S(a) [corresponding to H^°\x 7 ^) •■ 
(B8|) is satisfied for any function f(x) of x~ a (l — x) b type with a > 0, b > 0. For the model 

%{x 2 -a 2 



(38) 

c {x)\, the inequality 
y) which has a wider 



) profile and f(x) = A (1 — x) 3 , the inequality (38) is exactly saturated. If one takes the model (x, a) 



x f(x) S(x 2 — a 2 ) with an extremely wide profile, one obtains the result H^'(x,^) = ^{f(xi)/(l+£,) + f(x2)/(l— £)} 
which violates 



VI. SUMMARY 



In this paper, we treated double distributions as the basic objects for parametrizing nonforward matrix elements. 
An alternative description in terms of skewed distributions was obtained by an appropriate integration of relevant 
DDs. The use of DDs helps to establish important features of SPDs such as their nonanalyticity at the border points 
X = ( and x = ±£. DDs are crucial for securing the property that the moments of SPDs should be polynomial in 
the skewedness parameter. For these reasons, the use of DDs is unavoidable in constructing consistent models of 
SPDs. 
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